Key pairing interaction in layered doped ionic insulators 
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A controversial issue on whether the electron-phonon interaction (EPI) is crucial for high- 
temperature superconductivity or it is weak and inessential has remained one of the most challenging 
problems of contemporary condensed matter physics. We employ a continuum RPA approximation 
for the dielectric response function allowing for a selfconsistent semi-analytical evaluation of the 
EPI strength, electron-electron attractions, and the carrier mass renormalisation in layered high- 
temperature superconductors. We show that the Frohlich EPI with high-frequency optical phonons 
in doped ionic lattices is the key pairing interaction, which is beyond the BCS-Migdal-Eliashberg 
approximation in underdoped superconductors, and it remains a significant player in overdoped 
compounds. 
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For a long time, a basic question concerning the 
key pairing interaction in cuprate and other high- 
temperature superconductors has remained open. Some 
density functional (DFT) calculations [1, 2] found small 
EPI insufficient to explain high critical temperatures, Tc, 
while some other first-principles studies found large EPI 
in cuprates [s] and in recently discovered iron-based com- 
pounds [4]. It is a common place that DFT underes- 
timates the role of the Coulomb correlations, predict- 
ing an anisotropy of electron-response functions much 
smaller than that experimentally observed in the lay- 
ered high- Tc superconductors. The inclusion of a short- 
range repulsion (Hubbard U) via the LDA+U algorithm 
[5] and/or nonadiabatic effects [3] significantly enhances 
the EPI strength due to a poor screening of some par- 
ticular phonons. Also substantial isotope effects on the 
carrier mass and a number of other independent observa- 
tions (see e.g. [6] and references therein) unambiguously 
show that lattice vibrations play a significant although 
unconventional role in high-temperature superconduc- 
tors. Overall, it seems plausible that the true origin of 
high-temperature superconductivity could be found in a 
proper combination of strong electron-electron correla- 
tions with a significant EPI [7]. 

Here, we calculate the EPI strength with optical 
phonons, the phonon-induced electron-electron attrac- 
tion, and the carrier mass renormalisation in layered su- 
perconductors at different doping using a continuum ap- 
proximation for the renormalised carrier energy spectrum 
and the RPA dielectric response function. The Frohlich 
EPI with high-frequency optical phonons turns out the 
key pairing interaction in underdoped highly polarizable 
ionic lattices and remains significant player in overdoped 
compounds. 

We start with a parent insulator as La2Cu04, where 
the magnitude of the Frohlich EPI is unambiguously esti- 
mated using the static, and high-frequency, eoo dielec- 
tric constants To assess its strength, one can apply 
an expression for the polaron binding energy (polaronic 
level shift) Ep, which depends only on the measured 
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Here, the integration goes over the Brillouin zone (BZ), 
eo ~ 8.85 X 10"-^^ F/m is the vacuum permittivity, 
and hz = CsCQo/i^s — ^00) - In the parent insulator, the 
Frohlich interaction alone provides the binding energy of 
two holes, 2Ep^ an order of magnitude larger than any 
magnetic interaction {Ep = 0.647 eV in La2Cu04 [9]). 
Actually, Eq.([T]) underestimates the polaron binding en- 
ergy, since the deformation potential and/or molecular- 
type (e.g. Jahn- Teller [10]) EPIs are not included. 

It has been argued earlier [8] that the interaction with 
c-axis polarized phonons in cuprates remains strong also 
at finite doping due to a poor screening of high-frequency 
electric forces as confirmed in some pump- probe ml, ll2j 
and photoemission ^13.-.15j experiments. However, a 
quantitative analysis of the doping dependent EPI has 
remained elusive because the dynamic dielectric response 
function, e(cj,q) has been unknown. 

Recent observations of the quantum magnetic oscilla- 
tions in some underdoped [16] and overdoped [17] cuprate 
superconductors are opening up a possibility for a quan- 
titative assessment of EPI in these and related doped 
ionic lattices with the quasi two-dimensional (2D) car- 
rier energy spectrum. The oscillations revealed cylindri- 
cal Fermi surfaces, enhanced effective masses of carriers 
(ranging from 2me to Grrie) and the astonishingly low 
Fermi energy, Ep , which appears to be well below 40 me V 
in underdoped Y-Ba-Cu-0 [16] and less or about 400 
meV in heavily overdoped T12201 [17]. Such low Fermi 
energies [l8| make the Migdal-Eliashberg (ME) adiabatic 
approach to EPI [Tq] inapplicable in these compounds. 
Indeed, the ME non-crossing approximation breaks down 
at Xhw^/Ep > 1 when the crossing diagrams become 
important. The characteristic oxygen vibration energy 
is about huo = 80 meV in oxides [20, 21], so that the 
ME theory cannot be applied even for a weak EPI with 
the coupling constant A < 0.5. In the strong coupling 



regime, A > 0.5, the effective parameter XHujo/Ef be- 
comes large irrespective of the adiabatic ratio, hwo/Ep^ 
because the Fermi energy shrinks exponentially due to 
the polaron narrowing of the band [22]. Since carriers in 
cuprates are in the non-adiabatic (underdoped) or near- 
adiabatic (overdoped) regimes, Ep < huo^ their energy 
spectrum renormalised by EPI can be found with the 
familiar small-polaron canonical transformation at any 
coupling A 

The matrix element of the screened electron-phonon 
(Frohlich) interaction is found as [24] : 
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where Jo{q) is the bare (unscreened) vertex in the parent 
insulator. Fig. la. In our selfconsistent approach e(co', q) is 
calculated in the loop (RPA) approximation [25| with the 
exact (polaronic) carrier propagators taking into account 
the phonon "dressing" of carries: 
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Here, ek is the polaron band dispersion, r is the relax- 
ation time, Q is the normalization volume, and /k the 
Fermi-Dirac distribution function. The effect of collisions 
cannot be always taken into account merely by replac- 
ing by + z/r, in the collision- less Lindhard dielectric 
function, Eq.JS]), in particular at low frequencies or in the 
static limit [26, 27], where ladder-type vertex corrections 
are important [28j. However, in our case the relevant fre- 
quency is the (renormalised) optical phonon frequency, so 
that the vertex corrections are negligible as I/cjqt <C 1. 

Since the Fermi surfaces measured in the quantum os- 
cillation experiments [l6|, [l^ are almost perfect cylin- 
ders, one can apply the continuum (parabolic) approxi- 
mation for the quasi-2D polaron energy spectrum, ek = 
h'^{kl + ky)/2m* ^ where the polaron effective mass, m* 
has to be found self-consistently as a function of EPI. 
Calculating the sum in Eq.(|3]) yields the following di- 
electric response function extending the familiar pure 2D 
collision-less result by Stern [29^ to the 3D Coulomb in- 
teraction of carriers with the quasi-2D energy spectrum 
and collisions, 
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Here and below = -\- <?y is the square of the phonon 
momentum, N = kp/27TC is the carrier density, a and c 
are in-plane and c-axis (chemical) unit cell constants, re- 
spectively, and vf = hkF/m* is the Fermi velocity. The 
real and imaginary parts of susceptibility are found as 
(see also [30|) Xi(^,^||) = 2 - [R{z, z'^ - u, p) ^ R{z, z'^ ^ 



u,p)]/z^ and X2(^.^||) = [^(^ 
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FIG. 1: (Color online) Diagrammatic representation of the 
screened EPI vertex (a): solid lines correspond to polaron 
propagators, wavy lines are the exact phonon propagator, and 
the dashed line is the Coulomb repulsion, (b): long- wave dis- 
persion of zeros of the dielectric function of quasi- 2D carriers 
with the 3D Coulomb repulsion. An optical phonon with the 
energy 80 meV is also shown. 
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z = q\\/2kF^ u = uj/(2kFVF) and {3 = 1/(2/cf0 = ^ft 
is the mean- free path). 

As shown in Fig. lb, a collective mode, co'p^(q), de- 
fined by e(co'p^,q) = appears within the same fre- 
quency range as the optical phonon mode cjq, lead- 
ing to the plasmon-phonon mixing (so-called plasphons 
jsH). In the long wave limit we find using Eq.(j4]), 
^pi{^) ^ ^pq\\/q foi" collision- less carriers (/3 = 0), where 
hhjp = (e2£;F/47reoeooc)^/^ is c.a. 132 meV for Ef = 40 
meV . This mode is softer when it propagates across the 
planes than along the planes. Fig. lb, due to a low sus- 
ceptibility of quasi-2D carriers to the electric field applied 
across the planes [32]. 

The polaron level shift, Ep = V(0)/2, the carrier at- 
traction induced by EPI, — V(r), the in-plane polaron 
mass and the mass renormalisation exponent, are 
found as 
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FIG. 2: (Color online) The on-site, 1^(0) (upper curve), the 
nearest-neighbor, 1^(1) (middle curve), and the next-nearest- 
neighbor, V{V2) (lower curve), attractions induced by the 
Frohlich EPI for different doping x and two characteristic 
phonon frequencies, huoo = 80 meV and huoo = 60 meV. 
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respectively, where j connects nearest-neighbor sites. 
Holes in cuprates reside on oxygen, so that the nearest- 
neighbor hopping distances is j = a/>/2. The BCS 



FIG. 3: (Color online) The carrier mass (in units of the band 
mass m = rUe) as a function of doping for two characteristic 
phonon frequencies, huuo = 80 meV and huoo = 60 meV. 



coupling constant with phonons is defined as A = 
EpTTio? /jih? in the case of 2D carriers with a constant 
density of states {ma^ /2T:h? per spin), where m is the 
bare band mass in a rigid lattice. Using Ep > 0.6eV and 
m = 2me places cuprates in the strong-coupling regime, 
A > 0.86. 

Approximating the Brillouin zone as a cylinder of the 
volume 2t:'^ q^/c with the Debye momentum qo = 2y^/a 
and integrating over azimuthal (in-plane) angle yield 
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where = R{kt,t^/2 - u,P~p) + R{kt,t^/2 + u,P~p), 
^ = I{kt,t^/2 - u,k^) - I{kt,t^/2 + u.k^p), = 
e^c/(2T:'^eQKa^)^ O.TleV (with eoo = 5 and = 30 
of La2Cu04 [9|). The on-site, the nearest-neighbor 
and the next-nearest-neighbor attractions correspond to 
r = 0,l,v^, respectively. Jq{x) is the Bessel func- 
tion, k — kp/qo is the dimensionless Fermi momen- 
tum, u ujorn^o? /AttH ^ 0.015(m*/me) (for huj^ 80 
meV), C e^crrf /[eQKTT^h^) ^ 1.23(m*/me), and r] 
Ac^ /iio? ^3.93. If one assumes that carriers are scat- 
tered off impurities with the density equal to the carrier 
density, as in La2-a?Sra^Cu04, then in the Born approx- 
imation P = ^^{m* /rrie)'^ with Pq independent of the 
carrier mass and density. Using kpl ^ 20 and m* = 2me 
as found in YBasCuaOe.siie] yields /3q ^ 0.0125. 

We can now evaluate the doping {x) dependence of 
V{r) and m* by just changing the dimensionless Fermi 
momentum, k = (x/2)-^/^, in Eq.([TOj) and solving Eqs 
()10l8|9p selfconsistently with any bare band mass (we 
choose m = me in our numerical calculations). At very 
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low doping the on-site and nearest-neighbor inter-site at- 
tractions are enormous, V^(O) ~ 1.25 eV and V^(l) ~ 0.87 
eV, respectively, and carriers are rather heavy, m*/m ^ 
10. Such heavy polarons are readily localised by disorder 
accounting for the Mott variable range hopping, which 
explains conduction of lightly doped cuprates. 

With doping, the attraction and the polaron mass 
drop, Figs. (j2|3p . respectively. However, the on-site (r = 
0) and the inter-site (r = 1) attractions are well above 
the superexchange (magnetic) interaction J (about 100 
meV) in underdoped and optimally doped compounds 
since the non-adiabatic carriers cannot perfectly screen 
high-frequency electric fields. Both attractions and the 
mass renormalization remain also substantial at overdop- 
ing. The polaron mass, Fig.(|3|) agrees reasonably well 
with the experimental masses Increasing the 

phonon frequency enhances the attraction and lowers the 
polaron mass in underdoped compounds with little effect 
on both quantities at overdoping, Figs. (j2|3p . Decreasing 



(increasing) the band mass makes polarons lighter (heav- 
ier). 

In conclusion, we have quantified the carrier-carrier at- 
traction and mass renormalisation induced by EPI in lay- 
ered doped ionic lattices. The Frohlich EPI with high- 
frequency optical phonons turns out to be the key pairing 
interaction in underdoped cuprates and remains essen- 
tial at overdoping. What is more surprising is that EPI 
is clearly beyond the BCS-ME approximation since its 
magnitude is larger or comparable with the Fermi en- 
ergy and the carriers are in the non-adiabatic or near- 
adiabatic regimes depending on doping. Both conditions 
point to a crossover from the bipolaronic to polaronic 
superconductivity [22] with doping. 
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